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Abstract 

The production of dielectrons in reactions involving a\ mesons and pi- 
ons is studied. We compare results obtained with different phenomenological 
Lagrangians that have been used in connection with hadronic matter and fi- 
nite nuclei. We insist on the necessity for those interactions to satisfy known 
empirical properties of the strong interaction. Large off-shell effects in dielec- 
tron production are found and some consequences for the interpretation of 
heavy ion data are outlined. We also compare with results obtained using 
experimentally-extracted spectral functions. 



1 Introduction 

The field of heavy ion collisions is a very active one, straddling high energy and 
nuclear physics. At the upper energy limit of this flourishing area of research, the 
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goal is to eventually produce and study a new state of matter in the laboratory: 
the quark-gluon plasma (QGP). That strongly interacting matter in conditions of 
extreme energy densities undergoes a phase transition is in fact a prediction of QCD 
p]]. To confirm whether the phase transition indeed occurs in relativistic heavy ion 
collisions and that the QGP is formed, one needs a clear signal as a signature for 
the QGP. Many approaches have been suggested to elucidate the existence of this 
elusive state of matter, but unfortunately no single measurement can yet be singled- 
out as a "smoking gun" candidate. Instead, it appears that many complementary 
experimental data will require simultaneous analysis @. One class of observables 
that appears especially attractive is that of electromagnetic signals. This owes to the 
fact that such probes essentially suffer little or no final state interactions and thus 
constitute reliable carriers able to report on the local conditions at their emission 
site. Indeed, the calculated emission rates for photons and lepton pairs have been 
shown to strongly depend on the local density and temperature. Those facts were 
established some time ago and several experiments specializing in the measurement 
of electromagnetic radiation are either running right now or being planned. 

In this paper, we shall be concerned about lepton pairs and calculations of their 
production. With respect to photon measurements, they constitute an important 
related observable. The extra degree of freedom associated with the invariant mass 
makes them valuable tools in the study of annihilation reactions, for example. The 
emission of lepton pairs from the QGP and also from hot hadronic has been studied 
by many authors JJ. We will not discuss here the lepton pairs with high invariant 
masses (M > mj/y), they constitute a fascinating story of their own and the lit- 
erature on the subject is considerable. Rather, we shall address the issue of softer 
pairs and that of the reliability of calculations of their emission. Recently, a lot 
of attention has focused on measurements from the CERES/NA45 collaboration 
Q|. Several theoretical calculations have attempted to reproduce this experimen- 
tal data ||. Because of the exciting potential of such measurements to reveal new 
physics, extra care has to be taken in their theoretical modeling. In the very soft 
sector, relevant for the CERES experiment, different theoretical calculations for the 
fundamental processes have in fact been confronted wih each other with a striking 
degree of agreement ||. It has also been shown that this data supports the fact 
that the very soft sector of the lepton pair spectrum is relatively independent of the 
dynamical approach used to model the evolution of the four-volume of the colliding 
system. However, the disagreement between models and data has led to different 
interpretations |J. 

Another interesting aspect of the CERES analysis is that the region of invariant 
mass < M < mj/^ also shows an "anomalous" excess of lepton pairs, when 
compared with estimates based on hadronic electromagnetic decays. Note that this 
excess is seen only in nucleus-nucleus events, the proton-nucleus results are well 
reproduced assuming hadronic sources 0. Adopting the same line of reasoning 
as that used in most of the calculations done to interpret the low mass sector of 
the CERES experiment, it is natural to ask which hadronic reactions could con- 
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tribute significantly to the mass region in question. Some work on this aspect was 
performed in the past where it was found that in a hot environment of mesons, 
the reaction 7r + a\ — > e + e~ should be an important source of dileptons |7|, ||, || 
(from here on we shall discuss only dielectrons) . Such studies concerned with the 
contribution of the ai to the lepton pair spectrum are based on effective chiral La- 
grangians containing p, n and ai. Our purpose in this paper is the following: we 
review several different effective Lagrangians for the npai vertex, and we calculate 
lepton pair production from the 7r + a\ —>■ e + e~ reaction, using the Vector Meson 
Dominance model (VMD). We shall follow a philosophy similar to that in Ref. ||, 
but we go further by including a survey of other known and used interactions. We 
shall also stress the importance of conforming to empirical measurements, and we 
highlight the off-shell properties. Our work is organized as follows: after introducing 
the different models and mentioning their origins, we compare with each other the 
lepton results obtained and comment on the restrictions imposed on the effective 
interactions by the hadronic phenomenology. We also will show results from an 
experimentally-extracted spectral function analysis of lepton production processes. 
We then conclude. 



2 Effective Lagrangians and phenomenology of the 
a\ meson 

In our survey of the irpai interaction, we in turn will consider the following La- 
grangians: 

1. An effective chiral Lagrangian based on the SU(2)l x SU(2)r linear a model 

2. AU(2)lxU (2)r chiral model for pseudoscalar, vector, and axial- vector mesons 

3. An effective chiral Lagrangian where the vector mesons are introduced as mas- 
sive Yang-Mills fields of the chiral symmetry [|1^] . 

4. An effective Lagrangian previously used in connection with photon emission 
rates 



5. An effective Lagrangian used to address the issue of form factors in the Bonn 
potential |L4 . 



As one can realize, many interactions claim to address the a\^p vertex. This 
variety owes to the fact that there is currently no unique way of implementing the 
chiral symmetry in an approach uniting pseudoscalars, vectors and pseudovectors. 
Therefore, many versions of a given coupling can be obtained, with symmetry re- 
quirements sometimes yielding to specific practical concerns. Note that all of the 
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above Lagrangians (except one) have been used in the literature to calculate the 
lepton-producing reaction we are considering here. In such investigations, a pop- 
ular practice consists of adjusting any free parameter in the above interactions to 
reproduce the experimental decay width r ai In cases where more than one free 
parameter is available, one can also fit F ai ^^, when electromagnetic phenomena are 
of some relevance. However, another well-defined empirical quantity exists, regard- 
ing the decay ai —>■ np. The ratio of D-wave to S'-wave (eigenstates of the relative 
orbital angular momentum in the exit channel) amplitude has been measured and 
is known to be D/S = — 0.09 ± 0.03 ||15|| . Theoretically, the calculation of this 
ratio is accomplished by expanding the decay amplitude in the helicity basis, using 
the standard helicity representation for the polarization vectors. On the other hand, 
one can expand the amplitude in spherical harmonics. Working in the rest frame of 
the decaying particle, in this case the a%, one can then match the two representa- 
tions in terms of their helicity content. It is probably fair to say that this aspect of 
hadronic phenomenology has not received the attention it should have in the heavy 
ion community by users of effective hadronic Lagrangians. 

We will now proceed as advertised earlier and compare in turn the interactions 
described above in terms of their contribution to the process ir + a± — > e + e~ . How- 
ever, before we turn to this specific application, it is instructive to recall why the 
process in question is not determined uniquely by choosing an interaction and fitting 
on-shell properties. For our illustrative purpose, consider the process a\ — > pir. The 
most general vertex for this strong decay can be written as 

r"" = igi <T + i92 W + ig 3 Wq" + ig 4 + ig 5 WW , (1) 

where gi(k 2 , q 2 ) is a form factor, with W and q u being the the a% and p four-momenta. 
For on-shell a\ decay, three of the above terms are identically zero, owing to the 
transversality condition e(ai) • k = e(p) ■ q = 0. The on-shell form factors are 
usually chosen such that (some of) the experimental constraints discussed above are 
satisfied. However in the reaction under scrutiny here, n + a\ — > p — > e + e~ , the p 
meson is not on its mass shell and the extrapolations of form factors to this region 
are not unique. To evaluate the importance of such effects is the purpose of this 
work. Note however that off shell, for our reaction, there are still only two terms 
that are relevant in Eq. ([I]). This is because the lepton electromagnetic current 
is conserved. More specifically, the other term of Eq. ([I]) that could in principle 
contribute to our reaction is ig^fq v . However, the corresponding matrix element 
for electron-positron pair production, once the spin sums are carried out, will look 
like (the lepton masses have been set to zero, for simplicity) 

\M\ 2 cc gV^(7^i7^ 2 ) , (2) 

where q\{q%) is the four-momentum of the electron (positron). Note that now g M = 
Qi + 0.2 ■ The above expression, Eq. (H), is zero. 
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2.1 A Chiral SU(2) L x SU(2) R effective Lagrangian 

The effective Lagrangian for the a^p vertex constructed by Ko and Rudaz ]T0 



is 



based on the SU(2) L x SU(2) R linear a model, with the p vector meson and the 
ai axial-vector meson included as phenomenological gauge fields associated with a 
local SU(2)l x SU (2)r chiral symmetry. This local chiral symmetry is then broken 
to global chiral symmetry. 

The aipir vertex is generated by the following interaction Lagrangian |], |T0|j , 
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The ai(fc) — > 7r(p)p(g) decay width can be obtained as 
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where p is the momentum of the 7r meson evaluated in the rest frame of the a±, and 
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The ratio of D-wave to 5-wave amplitudes for this final state is 
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This chiral formulation admits a ain^y direct contact term and as such, devi- 
ates from the "traditional" Vector Meson Dominance approach. Including this, the 
radiative decay width one obtains for a\ — > 717 is 



r(ai7T7) 
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-p '-"ai 

Here is an opportune place for a short digression on the Vector Meson Dominance 
model (VMD). A recent discussion of its different representations can be found in 
Ref. |16| . We refer the reader to this reference, and those therein, for details. In 
short, the vector meson-photon vertex can either be a constant or the electromag- 
netic field strength tensor can couple to the rho field strength tensor to yield a 
g 2 -dependent vertex, where q is here the photon invariant mass [fT7|| . In the lat- 
ter case, an additional direct photon contact term is then necessary as the vector 
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Table 1 : Two parameter sets for the model of Ref . [|l0| , and the associated phenomenology. 



meson-photon mixing vanishes for real photons. The approaches described in this 
work invoke both versions of VMD. We have implemented the appropriate one for 
each case and also verified that gauge invariance was verified in the electromagnetic 
sector. 



Returning to the Lagrangian under scrutiny, following Refs. |L0], [18| we choose the 
parameters g,c,K&, and Z„ guided by phenomenological considerations. Two choices 
of those numbers are displayed in Table ^ Parameter set I was used previously to 
address dilepton production in the reaction we are concerned with here || . However, 
this choice produces a wrong sign for the D/S ratio and its \ 2 can be somewhat 
improved. This is achieved by parameter set II, see Table [l]. 

2.2 A U(2)lX U(2)p> effective chiral Lagrangian 

A £7(2) l x U(2)r chiral model for pseudoscalar, vector, and axial- vector mesons has 



been proposed by B. A. Li [11] which produces a successful description of meson phe- 
nomenology. In this effective chiral theory, the a\(k)n(p)p(q) coupling is described 
by the following interaction Lagrangian: 



C aip7T = Aa» ■ (p„ x t?) + BaT ■ {p v x d^ir), (8) 



where 
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In this model, f n = 0.186 GeV and the particle masses are taken as input. Fitting 
g = 0.35 yields a good list of light-meson empirical properties [|ll|].[] 



2.3 Another effective chiral Lagrangian 

Yet another effective chiral Lagrangian for pseudoscalar, vector, and axial-vector 



mesons can be derived [12]. In this model, the n meson is described through the 



nonlinear a model, and the p and a\ mesons are included as massive Yang-Mills fields 



of the chiral symmetry. This scheme has been used to discuss photon emission [p~9 
and dilepton production from hot hadronic matter. 

The lowest order interaction term for a\fm is given by a lengthy expression given 
in terms of the meson matrices |12[]. For the physical a± — > pir decay, the vertex 



function leading to the decay width is 



T^(ai — >■ p7r) = i(g ai p7rg^u - K^q^K) 
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Here, F n w 135 MeV. With the two sets of parameters which give the a± 



(14) 



pn 



decay width T(ai — > pn) w 400 MeV |TS| , we further obtained the D/S ratio for 
this effective Lagrangian: 



set I 
set II 
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10.3063, a 
6.4483, a = 



0.3405, £ = 
-0.2913, f 



0.4473, 
= 0.0585, 



D/S -- 
D/S 



0.357 
= -0.099 



(15) 



In this work, we shall use parameter set II. The electromagnetic interaction is in- 
troduced by imposing the U(l) em gauge symmetry on the effective chiral Lagrangian. 

1 In Eq.(62) of ref. [|Xl|] , we argue that a factor of y/2 is missing. For m ai — 1230 MeV, we then 
get D/S = -0.0895. The experimental value is -0.09 ± 0.03 ± 0.01. 



2.4 Two more simple effective Lagrangians for a\(m interac- 
tions 

Xiong, Shuryak and Brown |H| have defined an effective Lagrangian for aipn in- 
teractions, in order to calculate photon production from hot hadronic matter. This 
Lagrangian is 

£a lP TT = GpCL^ip ■ qg pu - q^Py)p v ii. (16) 

The coupling constant G p is determined by fitting the a\ —>■ pir decay width. Unfor- 
tunately, we can not reproduce the numerical results of Ref. fll3f . Using m ai = 1230 
MeV and fitting the total decay width, we get G p = 11.425 GeV -1 . Using the 
Vector Dominance Model we obtain the coupling constant of the aiTfj interaction 
as G 7 = 0.573 GeV^ 1 . This then implies a value of r ai _ )7r7 = 1.94 MeV, which is 
somewhat larger than the experimental measurement of 0.640 ± 0.246 MeV. This 
effective Lagrangian also predicts D/S = 0.185. 
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313.4 579.1 
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11.8 9.6 


fit 
0.067 
-0.099 
1.8 


331.7 
0.331 
-0.161 
3.1 


fit 
1.940 
0.185 
37.3 


fit 

0.312 
0.045 
7.3 


~ 400 MeV 
0.640 ± 0.246 MeV 
-0.09 ±0.03 



Table 2: A comparison of hadronic properties for the interactions discussed in this work. 
Note that the x 2 ' s appearing in this Table were calculated using the experimental mea- 
surements that appear on it only. When different parameter sets are involved for a given 
interaction, the distinction is explained in the text. 



The other simple effective Lagrangian which considers the aipn interaction was 
used to formulate a meson-exchange model for np scattering by Janssen, Holinde 
and SpethflTJ. It is 



C aifm = g ai (d,Pu - d v p„) ■ [tt x - (17) 

For m ai = 1230 MeV, the coupling constant g ai = 2.285 GeV -1 follows from the 
usual fit of the ai decay width. The D/S ratio here is D/S = 0.045. 

A comparative summary of the different on-shell properties and predictions is 
shown in Table 0. 



3 Results for dilepton production from hot hadronic 
matter 
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3.1 Effective Lagrangian approach 

In this section we shall compute rates of lepton pair emission with the effective 
interactions introduced earlier. In this study we rely on relativistic kinetic theory 
to provide an idealized dynamical framework. Because we are mainly interested in 
comparisons between different models, this line of reasoning is entirely adequate. 
We also set any chemical potential to zero, for simplicity. The temperature chosen 
for our baseline study is T = 150 MeV. 

In general, the dilepton production rate from the annihilation of two particles 
a(p\) and 6(^2) can be written asf2Ci[| 

d^dM 2 =U l ds i j2^j2^f fMfM dM 2 ^ i S - iPl+P2] 

(18) 

where M is an overall degeneracy factor, f(p) is the distribution function of the 
incoming particles at temperature T, v ie \ is the relative velocity of the two particles, 



J(Pi-P2) 2 -m 2 a m 2 b 

Vrel = =r-Z » ( 19 ) 

and a is the dilepton production cross section for the reaction a + b — > / + + Using 
Boltzmann distributions, one performs the phase space integrations and obtains 

f/A MtJttj Ki(M/T)X(s, ml , ml)a(M), (20) 



d 4 xdM 2 32tt 4 M iV ' ' v ' ai 
where K\ is a modified Bessel function, and A is the kinematic triangle function 

X(x,y,z) = x 2 -2x(y + z) + (y-z) 2 . (21) 

The cross section a(M) depends only on the square of the invariant mass of the 
dilepton s = (pi + P2) 2 = M 2 and is easily related to the square of a spin-averaged 
scattering amplitude. This in turn can be written as 



4f ) L^, (22) 



Where, a is the fine structure constant, and L^ u is a leptonic tensor given by 

L^ = q^ 2 +q^-g^ qi -q 2 , (23) 

and H^ u is a hadronic tensor for the process. Note that (qi + q^) 2 = s = M 2 . 
The hadronic tensor is calculated from the meson vertex function of the appropriate 
Feynman diagram with the relevant Lagrangian. 

We plot in Fig. [I] the rates for dielectron production from a gas of hot mesons, at 
a temperature of T = 150 MeV. We concentrate on irai reactions. All of the differ- 
ent interactions enumerated previously have been considered and their contribution 
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Figure 1: Dilepton production rates from the tt + a\ — > e + e reaction 
The curves labels are as follows: (a) Ref. Q, (b) Ref. @, (c) Ref. fll 
Ref. [10] parameter sets I and II (see Table pi), (e) Ref. [13). 



at T=150 MeV. 
(dl) and (d2) 



appears here. A striking feature is that the rates calculated with those Lagrangians 
span two-and-a-half orders of magnitude. Recall that all of them (except one) have 
been used in the literature to perform dilepton calculations very similar to the one 
done here. The features illustrated in Fig. |I| are essentially temperature- independent. 
We have verified this by also performing calculations at a temperature of T = 100 
MeV. We conclude from this that off-shell effects are indeed quite important. Con- 
sidering the curves labeled (dl) and (d2), one can verify by looking again at Table § 
that a modest change in the y 2 calculated on shell can result in an important vari- 
ation in off-shell behaviour. This tells us that although the on-shell \ 2 can be used 
as a goodness-of-fit parameter, it is far from being enough to specify unambiguously 
which interaction to use in situations like the one at hand. One has to proceed with 
caution. 

We now turn to an approach which should have the off-shell effects buit-in, in 
order to help us deciding upon an interaction to use. 

3.2 Using experiment ally- const rained spectral functions 

Above, we have found large differences in lepton producing rates using different 
hadronic Lagrangians to model the hard vertices. This result is worrisome, as many 
of those interactions have been used in the past for dilepton calculations, and there 
is no doubt that they will be used again in the future. We then turn to an alternative 
approach, with the hope that this will assist us in the selection of an appropriate 
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theory. In the context of heavy ion reactions, it has been argued recently that 
lepton pair production cross sections could in fact be determined from the spectral 
functions extracted from the inverse process: e + + e~ — > hadrons |2l|. We refer 



the reader to the appropriate reference for the complete details. It will suffice here 
to state that the rate for the emission of dilepton pairs of invariant mass M, at 
temperature T in a given reaction is given by 



dR 

dM* 



4rf 
~2tt 



MTK\{M/T) p cm (M) 



(24) 



where a is the usual fine structure constant, Ki(x) is a modified Bessel function, 
and p em (M) is a zero-temperature spectral function. This quantity is extracted from 
e + e~ annihilation data through 



p em ( S ) 



s a(e + e 



hadrons) 



16n 3 a 2 

We then follow the same approach as that of a current algebra calculation | 
extract the 7ia\ contribution from e + e~ — > An^ data. One can then obtain a(e H 
7rai), which can then be used in Eq. (|25|) to extract the spectral function. 
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Figure 2: The curve labels are the same as for the previous two figures. The double line 
represent the rates derived from cross sections calculated with the methods of Refs. [21, 22]. 



that this procedure is straightforward and relatively free from ambiguity, numerical 
or otherwise. It is tantamount to a direct evaluation of the dilepton cross sections 
by the same theoretical methods. 

We plot on Fig. the rates previously shown on Fig. [TJ, accompanied by the rate 
evaluated with the method outlined above. Note that this latter calculation does 
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not run over the entire invariant mass range covered in the plot, as we have chosen 
not to stretch the validity of the theoretical approach which was deemed optimal in 
the invariant mass range covered by the double line p2| . This coverage is sufficient 
for us to make our point. Considering that the spectral function is constrained 
by data, we use this analysis as an extra tool to discriminate between the various 
interactions considered in this work. We observe that the spectral function results 
are overshot by the interaction from Ref. 
rates related to the Lagrangians in Refs. 
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They are underpredicted by the 
and [TJJ], and by rates obtained with 
The interactions from Refs. 



parameter set II of the interaction from Ref. |10|| . The interactions from Refs. [T2 
and [|l(J (parameter set I) yield results in good agreement with the spectral function 
determination. Even though both Lagrangians seem to exhibit a satisfactory off- 
shell behaviour, going back to Table § one observes that the one from Ref. |12| also 
produces an excellent \ 2 f° r on-shell hadronic properties. It appears that, with the 
parameters used in this study and considering both on-shell and off-shell behaviours, 
this interaction achieves the compromise we were seeking. It is also a satisfying one 
from a theoretical point of view [p3[| . Note in closing that the rates above were 
determined in the narrow- width approximation for the a>\. While at first sight this 
limit seems unreasonable, it has been shown to affect little the magnitude of the 
thermal rates ||24||. Its main effect is to soften the threshold of this specific reaction. 



4 Conclusion 

We have calculated dilepton emission from nai reactions at finite temperature using 
several different Lagrangian found in the literature. We have found widely dif- 
ferent results. One might argue that some of those interactions were derived for 
entirely different purposes, therefore comparing them on the basis of dilepton emis- 
sion seems vaguely unappropriate. One should however remember that all of the 
Lagrangians (except that of Ref. |Tl|) have been used previously in such exercises. 
Up to now, a critical comparison of their results was lacking. With the help of an 
experimentally-constrained spectral function, combined with a quantitative analysis 
of on-shell properties, we were able to select an adequate interaction. It is clear that 
a companion study to this one will consider the rates for photon emission. This 
work is in progress and will be reported on later. There are however indications 
that for photons, the differences arising from the use of different Lagrangians will 
be less severe f8|. This probably owes to the fact that the mass shell condition for 
real photons has a restraining effect. 

It should be clear that our goal was not an exhaustive survey of the parameter 
space relevant to each of the models we have discussed. We viewed these interactions 
as representative of what is currently on the market. One could in principle devise a 
new completely phenomenological interaction. Our study highlights the important 
issue of what needs to be done in order to give credibility to results obtained in an 
environment as potentially complex as that of high energy heavy ion collisions. 
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